Model for Anisotropic Directed Percolation 



V. Lien Nguyen 
Institute of Physics, P.O. Box ^29, Bo Ho, Hanoi 1000 Vietnam 



Enrique Canessa* 

International Centre for Theoretical Physics, P.O. Box 586, 34100 Trieste, Italy 

(February 1, 2008) 



Abstract 



We propose a simulation model to study the properties of directed per- 
colation in two-dimensional (2D) anisotropic random media. The degree of 
anisotropy in the model is given by the ratio \x between the axes of a semi- 
ellipse enclosing the bonds that promote percolation in one direction. At 
percolation, this simple model shows that the average number of bonds per 
site in 2D is an invariant equal to 2.8 independently of \x. This result suggests 
that Sinai's theorem proposed originally for isotropic percolation is also valid 
for anisotropic directed percolation problems. The new invariant also yields a 
constant fractal dimension Df ~ 1.71 for all //, which is the same value found 
in isotropic directed percolation (i.e., [i = 1). 
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I. INTRODUCTION 



Critical phenomena Jl| in anisotropic systems without equivalent nearest neighbors con- 
stitute an interesting research topic [0 . A universal formula for percolation thresholds, that 
involves the dimension of the anisotropic lattice and an arithmetic average of the coordina- 
tion number for different anisotropic lattices, has been recently postulated in Ref. ||. The 
extension of these studies to more complex problems, such as directed percolation (DP), and 
more complex systems, such as anisotropic random systems, is yet to be addressed. In this 
context, random systems are good candidates to model anisotropy since they do not have 
equivalent nearest neighbors nor equivalent sites at all lengths. 

In this work we propose a simple simulation model to study the properties of DP in 
two-dimensional (2D) anisotropic random media. The degree of anisotropy is computed by 
means of the ratio fi = y c jx c between the axes of a semi-ellipse enclosing the bonds that 
promote percolation in one direction, such that y < x (see Fig.l). As a function of the order 
parameter fi and at the percolation threshold, we measure the correlation length exponent 
v and the fractal dimension Df of the largest percolating clusters (in systems of up to 51200 
random sites). In the present model, the well-known scaling exponents of isotropic DP follow 
by simply setting fj, — 1. 

At percolation threshold, our model shows that the average number of bonds per site for 
DP in anisotropic 2D random systems is an invariant (B c = 2.8) independently of fi. This 
result suggests that the Sinai theorem, proposed originally for isotropic percolation (IP), is 
also valid for anisotropic DP problems. The new invariant also yields a constant Df ~ 1.71 
for all fi, which corresponds to the value of isotropic DP. 

The paper is organized as follows. In the next section we outline our model. In Sec. Ill, 
we present the results of our simulations and discuss the effects of fi on the scaling exponents. 

II. THE MODEL 

In order to simulate DP in 2D anisotropic random media we develop a simulation al- 
gorithm similar to the one used in Ref. {§]. The coordinates of A sites are generated at 
random in a square box of size L = A 1 / 2 . The simulation length unit is chosen such that 
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the density of sites, namely n, in the box is always unity regardless of the total number of 
sites N. The percolation is then checked over sites from the left edge towards the right edge 
of the simulation box (i.e., along the x-axis in Fig.l). A periodical boundary condition is 
applied in the vertical y-direction. 

In Fig.l we show a 'particle' that moves from i to j. The moving is allowed whenever 
the site j is contained within the shaded elliptical area. In our simulations, the degree of 
anisotropy is given by the parameter \i = x m /y m , where x m is the longer and y m is the 
shorter axis of a semi-ellipse, i.e., [i is the ratio of the maximum 'hopping distances' along 
the x- and y-axes. 

In the standard 2D isotropic DP there are three possible equivalent directions to move: 
up, down and forward. This situation in our model is attained by setting u, — 1. In the limit 
fi — > oo, the model tends to the one-dimensional (ID) percolation problem. Thus, simulation 
results using the present 2D percolation model will reveal features of the crossover from the 
standard (say, isotropic) DP to the ID percolation problem. For intermediate values of 
fi > 1 our model features anisotropic DP. 

For a given value of the anisotropy parameter \i and for a given realization of random 
site coordinates {a^,^}, in a sample of size N = L x L, we study percolation from the left- 
to the right- simulation box edge. At the percolation threshold, we obtain the critical value 
of the semi-ellipse-axis y m : y c (N) = y m (N) and the mass M of the critical cluster: M(N) 
= " total number of sites belonging to the largest cluster at percolation" . These quantities, 
are then averaged over a great number K of random realizations of site coordinates for 
the same sample size N which result on the average quantities Y C (N) =< y c {N) > and 
M(N) =< M(N) >, respectively. 

In general, the dependence of the averages Y C (N) and M(N) on the samples size N is a 
consequence of the finite size effects of the percolation problem. In order to quantify these 
effects, the present simulations were performed at different N = 400, 800, 1600, 3200, 6400, 
12800, 25600 and 51200. Accordingly, the number K decreases from 10 4 to 10 2 such that 
the product of the numbers K x N is approximately the same for all sample sizes iV in our 
study. 

Along with these average quantities, we also calculate the moments 
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5y c (N) = < [y c (N) - Y C (N)] 2 >W 
5M{N) = < [M(N) - M(N)} 2 > 1/2 



(2.1) 
(2.2) 



and also the next-order moments, which are used to estimate the statistical errors of our 
simulation results. 

The present measurements are performed for various values of fi = 1,2,3,4,5,6,7 and 
8. As can be seen from the results discussed in the next section, the greater the value of /i, 
the stronger the finite size effects are. We verify that for \x = 8 simulations can only been 
carried out in samples of size N > 3200. 

Following the well-known finite-size scaling procedure suggested in Ref. ||, the critical 
exponent v of the percolation problem is defined from the scaling expression 

6y c (N) <x L- 1 /" , (2.3) 

where 5y c (N) is given in Eq. ( |2.1| ). Note that in the present study percolation is checked by 
the longitudinal direction only (the x-axes in Fig.l), then the exponent v in Eq.( |2.3| ) should 
be identified with the parallel v\\ (see |L|). 

III. RESULTS AND DISCUSSION 

In Fig. 2(a) the quantities — ki5y c (N) are plotted versus InL = lnA 1//2 for different 
values of the order parameter /i. The slopes of the fitting lines give the corresponding values 
for the exponent. Thus we found that for the largest \l — 8, v — 1.14 and for \x = 2 we 
measured v = 1.48. Other values are given in the figure. From these calculated moments 
and the linear fitting procedure, we estimate the statistical error to be less than 0.02 for all 
values of v shown in this figure. 

Results for the DP limiting case \l — 1 has been previously reported by one of us [|J. 
In this case, the value v « 1.65 ± 0.02 is known as the universal value of v\\ for a whole 
class of istropic DP models in 2D. As the amount of anisotropy increases, i.e. fi > 1, the 
correlation length exponent v decreases. Since this decrease is initially very fast to then 
become smoothly, it is not possible to obtain the whole crossover from 2D to ID directed 
percolation for the behaviour of v. That is, the decrease from v m 1.65 in 2D-DP to 
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the limit v = 1 in ID is limited by the size of our simulation box. The finite-size effects 
in correspondence to the different values of \i (and, therefore, to the different degrees of 
anisotropy in the 2D random systems) are in fact equivalent to those discussed in great 
detail in Refs. MM for anisotropic percolation and isotropic DP. 

By using the values of v in Fig.2(a), the critical 'radius' y c = y c {N — > oo) is determined 
from the scaling expression 

\y c -y c (N)\ cxL- 1 ^ . (3.1) 

In Fig. 3(a), (b) and (c) the quantities y c (N) are plotted versus L^ 1 ^ . From these plots we 
obtain y c by taking the asymptotic values iV — > oo for all fi studied. The estimated values 
of y c are also shown in this figure. 

Very remarkably, our simulations show that for all \i considered the quantity n x y 2 (fi) 
is in fact a constant. Since (ir/2)fiy 2 = (ir/2)x c y c is the area of the critical semi-ellipse at 
percolation, then our results suggests that Sinai's theorem proposed originally for IP, 
is also valid for 2D anisotropic DP problems. In this respect, we emphasize again that our 
length unit should be taken as n~ l l 2 for a system with site concentration n. 

Thus, our simulations leads to the new invariance 

Bf ] = S c n = (n/2)nfxy 2 = 2.82 ± 0.02 , (3.2) 

where n is the site concentration (e.g., the donor concentration in doped semiconductors), 
S c is the area of the critical semi-ellipse and is the mean number of connected bonds 



per site at percolation. The invariance of Eq^^) may be somehow related to the fractal 
behavior of the critical clusters as we shall discuss below. 

Let us determine first the fractal dimension Df of the critical percolation cluster using 
a standard procedure based on the scaling expression 

MiN)^ oc L°f. (3.3) 

In Fig. 2(b) the quantities In .M(iV) are plotted against InL for different values of the 
anisotropy parameter /x. 

Very surprisingly we found that the fractal dimensions Df, as determined from the 
slopes of the fitting lines for various values of fi in Fig.2(b), seem indeed to be constant and 



independent of /i within our simulation errors. We estimate T>f(fi) = Df pa 1.71 ± 0.02 for 
all /i, which corresponds to about the same value of the isotropic DP model with fj, — 1. 

At a first glance this result might rise some doubts, but we believe it can be understood 
in connection with the invariant given in Eq. ( |3.2| ). The invariance of B^ d \ with respect 
to changes in the anisotropy parameter [i, implies that the average number of connected 
bonds at percolation is independent of fi. If we assume the percolation process within 
an elementary semi-ellipse (as in Fig.l) to be the 'originating percolation rule', then the 
invariance of Eq.( [3.2|) could mean that the law to generate percolation clusters remains 
unchanged as /i varies. If this conjecture is right, we could suggest here a more general 
statement for all types of percolation models which are related to each other by the Sinai 
theorem; in these cases, the fractal dimensions of the percolation clusters could all be the 
same. 

It should be noted that our simulations are limited to 1 < \i < 8. It is in this range 
that we observed the invariance of B c and the constant value for Df. We believe that these 
features are maintained for a larger range of [i values. However, it is not feasible to increase \i 
and the sample size N simultaneously and get to the point where the present 2D simulation 
model crosses to the ID case (i.e., fi — > oo). 

To conclude, we have suggested a model for anisotropic directed percolation (ADP) and 
have presented the first simulation results for the main critical exponents of the model in 2D 
random systems. Quite surprisingly, we have found a new invariance for the average number 
of connected bonds at percolation due to presence of a suitable external force (e.g., shear 
stress, magnetic field, etc). Our simulations show that the product fi x y 2 c is a constant for 
all /i's considered. This invariance should be in close relation to the value of Df. 

We strongly believe the present model of ADP could be important to describe some 
physical phenomena such as hopping conduction in anisotropic n-Ge and n-Si under strong 
electrical fields, where the impurity wave functions are anisotropic and the conduction band 
splits into one ellipsoid [[|. Our measurements could be useful, for instance, in the ex- 
pressions for the hopping resistivity in 2D anisotropic random media. The new invariance 
~ 2.8 could be used in these systems similarly to the invariance B« = 4.5 for IP m a 
circle problem ||. We hope the present model will stimulate further investigations on this 
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direction. 
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FIGURES 



FIG. 1. The anisotropic directed percolation model. Percolation from site i to j is allowed 
whenever the site j is contained within the shaded elliptical area. The degree of anisotropy 
is given by the ratio /i = x m /y m , where x m is the longer and y m is the shorter axis of the 
semi-ellipse. The case fi = 1 yields the standard isotropic DP. At percolation, the average 
number of bonds per site (dots in this figure) is 2.8. 

FIG. 2. Simulation data and linear fitting for: (a) — \n5y c (N) of Eq. fl2.1| ) plotted against 
the sample size InL for several values of the anisotropy parameter /i; (b) \nAi(N) of 
Eq. ( |3.3| ) plotted against the sample size In L for several values of the anisotropy parameter 
fi. Within measurement errors, most fitted lines are parallel giving the average constant 
value Df ~ 1.71 ± 0.02. Note that fluctuations around this value are not systematic as fi 
increases. 

FIG. 3. Simulation data and linear fitting for y c (N) of Eq.( [3.1|) plotted against L~ x l v . 
(a) /i=l, 2 and 3; (b) /i=4, 5 and 6; (c) /i=7 and 8. The threshold y c for each fi studied is 
estimated from the asymptotic values iV — > oo. 
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